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Abstract — Relay networks having n source-to-destination pairs 
and m half-duplex relays, all operating in the same frequency 
band and in the presence of block fading, are analyzed. This 
setup has attracted significant attention, and several relaying 
protocols have been reported in the literature. However, most 
of the proposed solutions require either centrally coordinated 
scheduling or detailed channel state information (CSI) at the 
transmitter side. Here, an opportunistic relaying scheme is 
proposed that alleviates these limitations, without sacrificing the 
system throughput scaling in the regime of large n. The scheme 
entails a two-hop communication protocol, in which sources 
communicate with destinations only through half-duplex relays. 
All nodes operate in a completely distributed fashion, with no 
cooperation. The key idea is to schedule at each hop only a 
subset of nodes that can benefit from multiuser diversity. To 
select the source and destination nodes for each hop, CSI is 
required at receivers (relays for the first hop, and destination 
nodes for the second hop), and an index-valued CSI feedback at 
the transmitters. For the case when n is large and m is fixed, it 
is shown that the proposed scheme achieves a system throughput 
of m/2 bits/s/Hz. In contrast, the information-theoretic upper 
bound of (m/2) log log n bits/s/Hz is achievable only with 
more demanding CSI assumptions and cooperation between the 
relays. Furthermore, it is shown that, under the condition that 
the product of block duration and system bandwidth scales 
faster than log n log log n, the achievable throughput of the 
proposed scheme scales as (logn). Notably, this is proven to 
be the optimal throughput scaling even if centralized scheduling 
is allowed, thus proving the optimality of the proposed scheme 
in the scaling law sense. Simulation results indicate a rather 
fast convergence to the asymptotic limits with the system's size, 
demonstrating the practical importance of the scaling results. 

Index Terms — Ad hoc networks, channel state information 
(CSI), multiuser diversity, opportunistic communication, scaling 
law, throughput. 
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I. Introduction 

THE DEMAND for ever larger and more efficient wireless 
communication networks necessitates new network archi- 
tectures, such as ad hoc networks and relay networks. As such, 
there has been significant activity in the past decade toward 
understanding the fundamental system throughput limits of 
such architectures and developing communication schemes 
that seek to approach these limits. 

Among other notable recent results on the throughput scal- 
ing of wireless networks, Gowaikar et al. [1] proposed a new 
wireless ad hoc network model, whereby the strengths of 
the connections between nodes are drawn independently from 
a common distribution, and analyzed the system throughput 
under various fading distributions. Such a model is appropriate 
for environments with rich scattering but small physical size, 
so that the connections are governed by random fading instead 
of deterministic path loss attenuations. When the random 
channel strengths follow a Rayleigh fading model, the system 
throughput scales as logn. This result is achievable through 
a multihop scheme that requires central coordination of the 
routing of nodes. Moreover, full knowledge of the channel 
state information (CSI) of the entire network is needed to 
enable the central coordination. 

Along with the work on multihop schemes, such as [1] and 
[2], there is another line of work characterizing the system 
throughput for wireless networks operating with two-hop 
relaying. The listen- and-transmit protocol, studied by Dana 
and Hassibi [3] from the power-efficiency perspective, turns 
out to have interesting properties from the system throughput 
standpoint as well. This is in fact a two-hop amplify-and- 
forward scheme, where relays are allowed to adjust the phase 
and amplitude of the received signals. A throughput of 8(n) 
bits/s/Hz is achieved by allowing n source-to-destination (S- 
D) pairs to communicate, while m = 0(n 2 ) nodes in the 
network act as relays. It is assumed that each relay node 
has full knowledge of its local channels (backward channels 
from all source nodes, and forward channels to all destination 
nodes), so that the relays can perform distributed beamform- 
ing. Morgenshtern and Bolcskei worked in [4] with a similar 
distributed beamforming setup, and their results reveal trade- 
offs between the level of available channel state information 
and the system throughput. In particular, utilizing a scheme 
with relays partitioned into groups, and where relays in each 
group have CSI knowledge of only one backward and one 
forward channel, the number of relays required to support 
a O (n) throughput is m = Q(n 3 ). In other words, with 
lower level CSI, the number of required relays increases from 
0(n 2 ) to 0(n 3 ). An equivalent point of view is to state the 
throughput in terms of the total number of transmitting nodes 
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in the system, p = n + to. Then the system throughput is 
O (p 1 / 3 ), when the relays in each group know the channel 
for only one source-destination pair. When relays know the 
channels for all source and destination nodes, the throughput 
scales as (p 1 / 2 ). 

Although these works have made great strides toward under- 
standing wireless ad hoc network capacity, implementations of 
the schemes require either central coordination among nodes 
[1], [2] or some level of CSI (channel amplitude and/or phase) 
at the transmitter side [3], [4]. The centralized coordination 
between wireless relays does not come for free, since the 
overhead to set up the cooperation may drastically reduce the 
useful throughput [5]Q Likewise, in a large system, obtaining 
this level of CSI, especially at the transmitter side, may not 
be feasible. This paper addresses the need to alleviate these 
limitations by proposing an opportunistic relaying scheme that 
works in a completely decentralized fashion and imposes less 
stringent CSI requirements. 

A. Main Contributions and Related Work 

The main contributions of this work can be summarized as 
follows. 

• A two-hop opportunistic relaying scheme for operating 
over fading channels is proposed and analyzed. The 
scheme's salient features are: 

— It operates in a decentralized fashion. No cooperation 
among nodes is assumed or required. 

— Only modest CSI requirements are imposed. At each 
hop, each receiver is assumed to have knowledge 
of its local incoming channel realizations, while 
transmitters have access to only index-valued CSI 
via low-rate feedback from the receivers. 

• The throughput of the proposed scheme is characterized 
by: 

— It is shown that, in the regime of a large number 
of nodes n and fixed number of relays to, the 
proposed scheme achieves a system throughput of 
to/2 bits/s/Hz. This can be contrasted with the 
information-theoretic upper bound (to/2) log log n 
on the scaling of the throughput, achievable only 
with full cooperation among the relays and full CSI 
(backward and forward) at the relays. These results 
reveal an interesting feature of multiuser diversity: 
whereas full cooperation between relays can readily 
form parallel channels, and multiuser diversity can 
boost the throughput of each channel by a factor of 
log log n, when cooperation is not possible, multiuser 
diversity succeeds in restoring the parallel chan- 
nels, but must forsake the multiuser diversity factor 
log log n. 

— We show that m can grow (as a function of n) as 
fast as 0(logn), while still guaranteeing the linear 
throughput scaling in to. The linearity breaks down 

'However, in throughput scaling law studies, see, e.g., [1], [2] and [5], 
among many others, the overhead needed to set up cooperation is usually not 
explicitly accounted for. 



if m grows faster than O(logn). Furthermore, when 
the product of the block duration and the system 
bandwidth scales faster than log n log log n, the over- 
head due to feedback is negligible, and therefore, 
the achievable throughput scaling of the proposed 
opportunistic relaying scheme is given by 0(logn). 

• It is proven that, under the assumption of independent en- 
coding (i.e., no cooperative encoding) at the transmitters 
and independent decoding (i.e., no cooperative decoding) 
at the receivers, the system throughput is upper-bounded 
by the order of logn, even if centralized scheduling is 
allowed. This result is of interest in its own right, since 
it quantifies the system throughput of wireless ad hoc net- 
works under the scenario where neither transmitters nor 
the receivers can cooperate in avoiding and/or canceling 
interference. Thus, the network is interference-limited, 
unlike other works in which global CSI is assumed (and 
thus either cooperative encoding or decoding is possible), 
leading to a linear throughput scaling. The throughput 
scaling results under our pessimistic, yet more realistic, 
scenario, improve the understanding of throughput scal- 
ing of wireless ad hoc networks. 

• The proposed scheme is order-optimal in achieving the 
0(logn) throughput scaling. This suggests that, as far as 
throughput scaling is concerned, operating the network 
in a decentralized fashion, with local CSI at the receivers 
and low-rate feedback, incurs no loss. 

• Simulation results show that the asymptotic conclusions 
developed in this work settle rapidly. Hence, the above 
scaling laws provide rule-of-thumb guidance for the de- 
sign of practical wireless systems. 

The key idea behind the proposed scheme is to schedule 
at each hop only the subset of nodes that can benefit from 
multiuser diversity. The concept of multiuser diversity was 
originally studied in the context of cellular systems [6]-[8]. It 
is known that the capacity of single-cell system is maximized 
by allowing only the user with the best channel to transmit 
at any given time. The concept is by now well understood 
in the context of infrastructure wireless networks, and has 
been adopted in 3G cellular systems and other emerging 
wireless standards [9]. However, to the best of our knowledge, 
it has received less attention for wireless ad hoc networks, 
with some exceptions such as [10], in which the potential 
of opportunistic relaying is reported in a setup with one 
S-D pair and multiple relay nodes, and the focus is on 
diversity-multiplexing trade-off analysis [11]. In this work, we 
highlight another aspect of multiuser diversity: its application 
to simplify network operations and its effect on throughput 
scaling. The opportunistic scheme proposed here is in the spirit 
of [12], where distance-dependent, random channel gains were 
exploited in scheduling. 

In this work, we restrict ourselves to those assumptions 
that are implementable with the state-of-the-art technologies. 
Specifically, we focus on the assumptions of perfect CSI at 
the receivers and partial CSI at the transmitters via low-rate 
feedback. With these less idealistic CSI assumptions, it is 
envisioned that independent encoding at the transmitters and 
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independent decoding at the receivers are employed. To the 
best of the authors' knowledge, there are no analogous results 
in the literature that consider the same scenario. It was recently 
shown, however, that under more optimistic assumptions on 
CSI in the network, linear throughput can be achieved by 
either joint encoding at the transmitters [13] or hierarchical 
cooperation with joint decoding at the receivers [5]. 

B. Organization of the Paper 

The rest of the paper is organized as follows. The sys- 
tem model and the proposed two-phase relay protocol are 
introduced in Section [TT] Section [III] characterizes the sys- 
tem throughput in the regime of large n and fixed m. The 
throughput scaling of the proposed scheme is evaluated in 
Section [IV] by explicitly taking the feedback overhead into 
account. Also in Section IIV1 the throughput upper bound is 
developed, valid even when centralized scheduling is allowed. 
Section [V] presents numerical performance results. Section [VTl 
briefly discusses the impact of relay cooperation on system 
throughput and the delay consideration. Finally, Section IVIII 
concludes the paper. Technical details and proofs are placed 
in the appendices. 

Notation: The symbol \X\ denotes the cardinality of the 
set X, and unless specified otherwise, log(-) indicates the 
natural logarithm. We write X ~ Exp(l) to indicate that 
the random variable X follows the standard exponential 
distribution with probability density function (pdf) given by 
fx(x) = exp(— x), x > 0. The indicator function is denoted 
by l(-), and we use "x 2 (2p)" to denote a chi-square random 
variable with 2p degrees of freedom. For two functions f(n) 
and g(n), f(n) = 0(g(n)) means that lim n _,oo|/(n)/g(n)| < 
oo, and f(n) = Sl(g(ri)) means that g(n) = 0(f(n)). We 
write f(n) = o(g(n)) to denote lim n _ KX3 |/(n)/5(n)| = 0, 
and f(n) = Q(g(n)) to denote f(n) = 0(g(n)) and 
g(n) = 0(/(n)). * 

II. System Model 

Consider a wireless network with n S-D pairs and m relay 
nodes, all operating in the same frequency band of width W 
Hz, in the presence of fading. Ad hoc nodes that generate data 
traffic are referred to as source nodes; nodes that receive data 
traffic are referred to as destination nodes. Relay nodes have 
no intrinsic traffic demands. We consider a two-hop, decode- 
and-forward communication protocol, in which sources can 
communicate with their destinations only through half-duplex 
relays. In the first hop of the protocol, a subset of sources 
is scheduled for transmission to relays. The relays decode 
and buffer the received packets. During the second hop of the 
protocol, the relays forward packets to a subset of destinations 
(not necessarily the set of destinations associated with the 
source set in the first hop). These two phases are interleaved: 
the first hop is run in the even-indexed time-slots, and the 
second hop is run in the odd-indexed time-slots. An example 
of the two-hop relay protocol is depicted in Fig. [T] 

We first describe the channel model. It is assumed that the 
wireless network has independent and identically distributed 
(i.i.d.) Rayleigh connections hi, r from source nodes i, 1 < i < 



n, to relay nodes r, 1 < r < m. Thus, the channel gains follow 
an exponential distribution, i.e., 7;. r = |/i; ir | 2 ~ Exp(l). 
Likewise, we assume that the channel gains £ rj from relays 
r, 1 < r < m, to destination nodes j, 1 < j < n, 
be i.i.d. Exp(l), and that channel gains ji, r and £ r j are 
independent for all i, r, and j. This model is appropriate 
for dense networks in a rich scattering environment, where 
the distance between transmitters and receivers has only a 
marginal effect on attenuation, and the channel attenuation is 
dominated by the small-scale fading due to multipath. Quasi- 
static fading is assumed, with channel gains fixed during 
the transmission of each hop, which is assumed to have a 
duration of T seconds, and taking on independent values at 
different transmission times. In practice, T can be as large 
as the coherence time of the channel allows. Regarding CSI, 
we assume that at each hop, each receiver has perfect CSI 
knowledge, while the transmitters have access only to an index 
value via receiver feedback used to indicate a source chosen 
for transmission. This CSI assumption is reasonable in practice 
as most wireless access network standards incorporate some 
form of pilot signals, and the type of feedback specified has 
low overhead. 

We now describe the scheduling at each hop. We start 
with the first hop (Phase 1). All relays operate independently. 
Thus, without loss of generality, let us focus on any specific 
relay, say r. By assumption, relay r has the knowledge of 
7i r , i = 1, . . . , n, and it will schedule the transmission of the 
strongest source node, say i r = argmaxj 7;, r , by feeding back 
the index i r at the beginning of the block. The overhead of 
this phase of the protocol is a single integer per relay node. 
Suppose the scheduled nodes constitute a set JC C {1, . . . , n}; 
then since there are m relays, up to m source nodes can be 
scheduled in this fashion, i.e., \JC\ < m. It is noted that it is 
possible for multiple relays to schedule the same source. In 
such cases, \JC\ < in. The scheduled source nodes transmit 
simultaneously with constant power P and fixed transmission 
rate of 1 bit/s/Hz@ Each relay sees a superposition of all the 
transmitting signals, i.e., 

y r = \fPhi T , r Xi r + V~P h t , r x t + n r , r = 1, . . . , m 
teic 

where Xi denotes the transmitted signal of source i, n r denotes 
the additive noise at relay r. In this paper, we assume Xi's 
are letters from codewords of a Gaussian capacity-achieving 
codebook satisfying E[|xi| 2 ] = 1. We further assume n r 's 
are i.i.d. complex Gaussian with zero mean and unit variance 
CAf(0, 1) and are independent of the fading channels. Since 
the transmission rate is 1 bit/s/Hz, communication can be 
supported in an information-theoretic sense if the correspond- 
ing signal-to-interference-plus-noise ratio (SINR) is greater or 
equal to one, i.e., 

SINR^ r = — — ^ >1, r=l,...,m (1) 

2 Generalizing to higher transmission rates is straightforward, but it encum- 
bers notation without adding insight. See discussions in Remark [6] for the 
motivation of choosing the rate of 1 bit/s/Hz. 
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Fig. 1. A two-hop network with n = 5 S-D pairs and m = 3 relay nodes (denoted by the blue disks), (a) In the first hop, source nodes {3, 4, 5} transmit 
to the relays, (b) In the second hop, the relays transmit to the destination nodes {2,3,4}. Solid lines indicate scheduled links, while dashed lines indicate 
interfering links. 



where p = P is the average signal-to-noise ratio (SNR) of the 
source-relay link. 

The scheduling at the second hop (Phase 2) works as 
follows. All relay nodes transmit simultaneously with fixed 
power Pr. Assume that the additive noise at the destination 
nodes are i.i.d. CAf(0, 1) and are independent of the fading 
processes {£, r .j}, and assume that independent messages are 
sent across relay nodes (which is the case in the proposed 
scheduling), destination node j can compute m SINRs by 
assuming that relay r is the desired sender and the other relays 
are interference as follows: 



SINR^j = 



1/ PR + El<<<m &,j ' 



l,...,m (2) 



where p^ = Pr denotes the average SNR of a relay- 
destination link. If the destination node j captures one good 
SINR, say, SINR^ > 1 for some k, it instructs relay k to 
send data by feeding back the relay index k at the beginning 
of the block. Otherwise, the node j does not provide feedback. 
It follows that the overhead of the second hop is also at most 
an index value per destination node. When scheduled by a 
feedback message, relay k relays the data to the destination 
node at rate 1 bit/s/Hz. In case a relay receives multiple 
feedback messages, it randomly chooses one destination for 
transmission. It is noted that in steady-state operation of the 
system, the relays have the ability to buffer the data received 
from source nodes, such that it is available when the oppor- 
tunity arises to transmit it to the intended destination nodes 
over the second hop of the protocol. This ensures that relays 
always have packets destined to the nodes that are scheduled. 
In addition, due to the opportunistic nature of scheduling, the 
received packets at the destinations are possibly out of order, 
and therefore each destination is assumed to be able to buffer 
them before decoding. 

Remark 1: It is noteworthy to draw a comparison between 
Phase 1 and Phase 2. From the relays' perspective, both hops 
of the communication protocol rely on scheduling a subset of 
"good" source/destination nodes for transmission. However, 
these two phases of the protocol differ in one key aspect: 
transmission over the second hop can be guaranteed to be 
successful since receivers (the destinations) have access to 



the SINRs, but this is not the case for the first hop. This is 
because in the first hop, each receiver (relay) selects a source 
node without knowledge of what the other relays select. As 
a consequence, each relay has no access to the interference 
stemming from all other concurrent transmitting sources, and 
therefore has no a priori knowledge of its own SINR. For 
example, in (Q~|i relay r knows the desired link strength 7< r , 
but it does not know JC and the corresponding interference 

term 2~2teKt^i 7*>f F° r trie seconc l no P> the senders (now 
the relays) are known a priori, and therefore the destination 
nodes have direct access to SINRs. This implies that once the 
destination node captures an SINR P2 > 1, and accordingly 
requests a transmission, this transmission will be successful at 
a data rate of 1 bit/s/Hz. This key difference between the two 
phases is mirrored in the analysis in Section [HI] 

Remark 2: In both hops, independent encoding at the trans- 
mitters and independent decoding at the receivers are em- 
ployed. By independent encoding, it is meant that the transmit- 
ters encode their message independently. This is a consequence 
of the fact that the transmitters have access to only partial CSI. 
Similarly, by independent decoding is meant that receivers 
decode their message independently by treating interference as 
noise. It is worth pointing out that, with the assumption of CSI 
at the receivers, techniques like interference cancelation are 
possible at the receivers. However, as elaborated in Remark|6] 
they are not interesting in our setup and thus are not considered 
here. 



III. Throughput: Large n and Fixed m 

Motivated by the observation that as communication devices 
(source and destination nodes in our system) become more and 
more pervasive, the number of infrastructure nodes (relays) is 
not likely to keep pace, the throughput analysis in this paper 
pays special attention to a regime in which the number of 
source and destination nodes, n, is large, while the number 
of relay nodes, to, is relatively small. We show that both 
Phase 1 and Phase 2 achieve average throughput (by averaging 
over random channel gains) of m bits/s/Hz, yielding am/2 
bits/s/Hz throughput for the complete two-hop scheme. We 
also show that for any two-hop protocol, the throughput is 
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upper-bounded by ^ log log n bits/s/Hz. This information- 
theoretic upper bound holds even if we allow full cooperation 
between relays and assume full CSI is available at the relays. 
Thus, the proposed scheme, with much simplified assumptions 
of decentralized relay operations and CSI at the receiver, 
succeeds in maintaining the linearity of the throughput in the 
number of relay nodes. 

A. Phase 1: Source Nodes to Relays 

In Phase 1, to relays operate independently and each 
schedules one source node for transmission. Hence, the total 
number of scheduled source nodes can be any integer between 
1 and to, i.e., \JC\ < to. In cases when \IC\ < to, multiple 
relays schedule the same source node, and the analysis of 
the probability of successful transmission should consider 
explicitly those links with multiple receivers. Due to the mul- 
tiplicity of possible combinations, the exact characterization 
of the average throughput of Phase 1, Ri, is mathematically 
involved. Fortunately, in order to show the achievability of m 
successful concurrent transmissions, it suffices to lower-bound 
i?i by considering only cases in which the m scheduled source 
nodes are distinct (thereby discarding the contributions to the 
throughput of the other combinations). 

By symmetry, each source node has a probability of 1/n to 
be the best node with respect to a relay. Thus, the probability 
that the scheduled users are distinct, i.e., no source node is 
scheduled by more than one relay, is given by Pr[7V m ] = 
n(n — 1) ■ • ■ (n — m + l)/n m , where N m denotes the event 
{to distinct source nodes are scheduled}. Now, a lower bound 
on R\ is 

m 



R 1 > Pr[N m ] ^Pr[SINR ? „, r > l] • 1 

r=l 



to Pr[iV m ] Pr[SINR pl > l] 



(3) 



where, for notational brevity and by the i.i.d. channel model, 
we drop the source node and relay indices in the last equation 
and use SINR P1 to denote the SINRs at all relay nodes. 

Now, we focus on the Pr[SINR pl > l] term. Again, for 
notational convenience, for a realization of n i.i.d. random 
variables X\, . . . , X n , we introduce X := m&x{Xi, . . . , X n } 
and Y := YlieK. 1 ^ i wnere is anv randomly selected (to — 
l)-element subset of {1, . . . ,n}\{j : Xj = X}. With these 
definitions, we have Pr[SINR pl > l] = ^[j/ffy ^ For 
the Rayleigh fading case, in which the link strengths are i.i.d. 
Exp(l) random variables, the cumulative distribution function 
(cdf) of X (largest of n i.i.d. Exp(l) random variables) can 
be written explicitly as Fx(x) = (1 — e~ x ) n . The asymptotic 
properties of X are well studied in literature (see [8] and [14]). 
For so = log n — log log n, it can be shown that Pr[X < sq] — > 
[14, eq. (A4)]. With the help of this property, we proceed to 
lower-bound Pr[SINR pl > l] by introducing a real variable 
s (0 < s < sq). By the law of total probability, we have 



Pr 



X 



1/p + Y 



> 1 



= Pr[X > s] ■ Pr 
+ Pr[X < s] • Pr 



X 



1/p + Y 
X 



> 1 



X > s 



1/p + Y 



> 1 



X < s 



> Pr[X > 

> Pr[X > 
= Pi[X > 



Pr 



Pr 



Pr 



X 



1/p + Y 



1/p + Y 



> 1 



> 1 



X > 



X > 



1 



1/p + Y 
F x (s))Fy(s-l/p), 



> 1 



where (HJi follows from the fact that, 
Pr[X < s ] 0, we have Pr[X 
p r [_ T JL^ > \\x > si -> Pr 



with 

> *] 
> 



< 



1] 



(4) 

(5) 

s < so and 
1, and thus 

.1^ ^ -I- - -J - x 1 iTJp+Y ^ -J ( which can 
be trivially shown by the law of total probability). Note that 

the lower bound (O suggests a suboptimal scheduling scheme 
according to which, each relay schedules the transmission of 
the "strongest" source only if the source's power gain exceeds 
a prescribed threshold s. The probability of such event is given 
by 1 — Fx(s), and Fy(s — 1/p) is a lower bound on the 
probability of a successful communication with the relay at a 
rate of 1 bit/s/Hz. 

The characterization of distribution of the interference term 
Y in © needs more care. This is due to the fact that, condi- 
tioned on not being the maximum among n channel strengths, 
each interference term in Y is no longer exponentially dis- 
tributed, and the interference terms are not independent in 
general. However, as shown in Appendix |A] these properties 
hold asymptotically with n. Numerical results in Appendix [A] 
show that these asymptotic trends are achieved for relatively 
small values of n. Thus, we can approximate Y as chi-square 
random variable with 2(m — 1) degrees of freedom, whose cdf 
is thus given by [15, eq. (2.1-114)] 



F Y (y) = l-e-yJ2 



k=0 



y 



(6) 



Substituting © and © into (O yields the following lower 
bound on the throughput of Phase 1. 

Lemma 1: For any p, m and < s < log n — log log n, 
the achievable throughput of the opportunistic relay scheme 
in Phase 1 is lower-bounded by 

, n(n-l)---(n~m+l) _ ^ — 



Rl > TO- 



1 



-(s-l/p) 



m-2 

E 

fc=0 




(7) 



as n — * 00. 

A tighter lower bound can be found by maximizing (f5]) 
over s, but we find that little insight can be gained from this 
exercise. The tightness of the lower bound (O is substantiated 
by numerical results shown in Fig. [2] of Section [VI 

Remark 3: Inspecting (0, we note that the lower bound on 
Ri exhibits a tradeoff between quantity and quality of sched- 
uled links. By increasing the number of relays to, one can 
schedule more simultaneous transmissions, which is beneficial 
from the throughput perspective. However, more transmissions 
generate more interference, degrading the SINR and lowering 
the probability of successful transmissions. In fact, as shown 
in Section |Vj not only the lower bound discussed here, but 
also the actual throughput i?i demonstrates this tradeoff. The 
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characterization of the best to (in terms of scaling) that 
maximizes throughput is pursued in Section [IV] 

For the regime of interest, where n is large and m fixed, it 
can be trivially shown, e.g., by setting s = logn — log log n, 
that the above lower bound approaches to. Note that this is 
also the best we can hope for in Phase 1, since the SINR > 1 
constraint to decode a transmitter implies that no more than 
to sources can be successful. 

The following corollary to Lemma Q] follows immediately. 

Corollary 1: For fixed to, Ri — > to as n — > oo. 

B. Phase 2: Relays to Destination Nodes 

We now develop an exact expression for the sum-rate of 
the relay-destination links. This is done by first showing that 
only a single relay per destination can produce a required 
SINR larger than one, and then computing the probability of 
the event that a relay is scheduled and consequently delivers 
throughput. 

Lemma 2: For any pn, to and n, the achievable throughput 
of the opportunistic relay scheme at Phase 2 is given by 



R 2 = m 1 - 1 - 



-l/PK 



(8) 



Before embarking on the proof, it is worth examining the 
statistics of the SINR in (f2]i. Given the i.i.d. channel model 
introduced in the previous section, the SINRs measured at 
each destination (cf. (f2]i) are of the generic form SINR^ = 

, , ftfl With the help of ©, the pdf of the SINR can 

l/PH+X (2m-2) r 1—1" r 

be shown as [14]: 



f(x) = / f(x\y)f Y (y)dy 
Jo 

e -X/p R / I 



(l + s)*' \PR 

The corresponding cdf is 



(1 + x) + m — 1 



F{x) = 1 - 



-x/PR 



1 + x) 



x>0. 



(9) 



(10) 



Note that the SINR r -s are i.i.d. over j = 1, . . . , n (but are not 
independent over r = 1, . . . , rri). 

Proof: First, we observe that each destination node j has 
at most one SINR^ > 1 for all relays 1 < k < to. To see this, 

P2 

assume SINR^ > 1 for some relay k, and consider another 
index k' 7^ k. From (O, we have 

l<£<rn 

from which it follows that 



Therefore 



SINRj^j = 



k',j 



I/PR- 



2^2l<t<m £,1,3 



< 



6=j 



< 1. 



Thus, each destination node can have at most one good relay 
as its sender. 

Now the sum-rate for Phase 2 depends on how many relays 
are scheduled by destinations. The probability that a relay finds 
no destination satisfying SINR > 1 i^| 

Pr[relay r does not receive feedback] 



Pr 



SINRPj < 1, Vj 



r,3 



(F(1)T 



-1/PR 



1 



The throughput of the relay-destination links is given by 
summing the probabilities of the relays engaged in transmis- 
sion. Accounting for the 1 bit/s/Hz rate per relay, we have that 
the average throughput of the second hop is given by 

m 

i? 2 = Pr [relay r transmits data to a destination] ■ 1 




(F(i)y 



1-1- 



,-l/PR 



(ID 

I ■ (12) 

This completes the proof of the lemma. □ 
The following corollary ensues by direct computation: 

Corollary 2: For fixed m, R2 — > to as n — » 00. 

Remark 4: It is interesting at this point to draw a connection 
between the scheduling of Phase 2 of the opportunistic scheme 
proposed here with the random beamforming scheme due to 
Sharif and Hassibi in the context of multiple-input multiple- 
output broadcast channels (MIMO-BCs) [14]. Seemingly un- 
related, the SINRs of both setups turn out to have the same 
distribution (cf. (O). To explain this subtlety, note that in the 
random beamforming scheme of [14], a random unitary matrix 
is applied to the data streams before sending them over the 
channel H (hence the terminology "random beamforming"). 
With the assumption of i.i.d. Rayleigh fading, entries of H 
follow i.i.d. circularly symmetric complex Gaussian random 
variables CAf(0, 1). By the isotropic property of the i.i.d. 
complex Gaussian random matrix H, #ff has the same 
distribution as H [16]. It follows that the channel statistics 
of the random beamforming scheme in the beam domain are 
the same as in the original antenna domain. In other words, 
the SINR in the beam domain is still of the generic form 



SINR = 



X 2 (2) 



(13) 



l/p fl + X 2 (2TO-2)' 

which is the same as in our Phase 2 (cf. (f2]i). 

Despite the mathematical equivalence, our proposed scheme 
for Phase 2 simplifies the random beamforming scheme in 
several respects: 

3 It may seem logical to turn off a relay for which the highest SINR is 
still less than one, but we still allow such relays to transmit (say, control 
information). This is because, as shown by numerical results, the performance 
is limited by the source-relay link. 
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• Random beamforming requires cooperation among the 
transmitters to form a beam. Opportunistic relaying op- 
erates in a completely decentralized fashion. 

• Random beamforming requires the feedback of an integer 
(the beam index) as well as a real number (the in- 
stantaneous SINR). The proposed opportunistic relaying 
scheme requires the feedback of only an index number. 
This simplification is justified by [14, Th. 2], which 
implies that when the system operates in the limit as 
n — ► oo with m = O(logrt), the aggregate interference 
from concurrent transmissions eventually hardens the 
instantaneous SINR near the value 1. Thus, there is no 
longer a need to feed back the SINR value. Furthermore, 
in terms of the throughput scaling law (as discussed later 
in Section IPVT ). this simplification incurs no loss. 

C. Feedback Overhead Analysis 

A detailed study of feedback overhead in the regime of 
large n and fixed to is omitted here for the sake of brevity. 
The calculation can follow the same steps as in Section HV-CI 
where we present a detailed analysis of the feedback overhead 
in the limiting regime of large n and m. 

D. Two-Hop Communication 

With the help of Corollaries Q] and |2j and by taking into 
account the 1 /2 penalty due to the two hops, the overall system 
throughput, defined as | min{-Ri, R^}, can be readily shown 
to be given as follows. 

Theorem 1: For fixed to, the two-hop opportunistic relaying 
scheme achieves a system throughput of m/2 bits/s/Hz asn-> 
oo. 

Since the proposed scheme works in a decentralized fashion 
and with low rate CSI feedback, it is natural to expect some 
throughput degradation compared to more intensive schemes. 
We will show that the opportunistic relaying scheme exhibits 
the pre-loglog factor of the scaling law of the throughput of 
more intensive schemes. To see this, we find an information- 
theoretic upper bound on the achievable scaling law for the 
aggregate throughput of any two-hop relaying scheme. 

Lemma 3: For any two-hop relaying architecture, with 
fixed to and SNR, the sum rate capacity scales at most as 
imloglogn as n — > oo. 

Proof: In two-hop relay schemes, all data traffic passes 
through relays. Therefore, the best scheme would be one in 
which all m relay nodes can cooperate and the relays have full 
CSI (i.e., backward as well as forward channel realizations). 
In such case, the two-hop communication can be interpreted 
as MIMO multiple access channels (MACs) followed by 
a MIMO-BC. The capacity region of the MIMO-BC, and 
the optimality of dirty-paper-coding (DPC) in achieving the 
capacity region have been shown in [17]. Furthermore, the 
capacity scaling of the DPC scheme is shown in [14] to be 
m log log n, which is also the capacity scaling for MIMO- 
MAC due to the MAC-BC duality [18]. Now, Lemma [3] 
follows by taking the two-hop penalty 1/2 into account. □ 

Remark 5: Contrasting Theorem [T] to Lemma [3] reveals 
two different facets of multiuser diversity. Fundamentally, 



multiuser diversity gain is a power gain, e.g., in the Rayleigh 
fading case, multiuser diversity schedules the best user for 
transmission, and boosts the average power by a factor of log n 
[8]. With the assumption of relay cooperation, as in Lemma[3] 
a spatial multiplexing gain equal to the number of relays to can 
be readily achieved (e.g., even by a suboptimal zero-forcing 
receiver [16]). Then, multiuser diversity can further boost 
the rate of each parallel channel by log log n, as shown by 
Lemma[3] In contrast, with the proposed opportunistic scheme, 
where relays operate independently, there is no guarantee 
of achieving the multiple parallel channels. Here, multiuser 
diversity is used as a mechanism that compensates for the 
interference plus noise so that the scheduled link can support 1 
bit/s/Hz. Ultimately, one achieves the linear scaling in to. Note 
that only with multiuser diversity gain does the SINR of each 
noncooperative link have the chance to meet the threshold^ 

Remark 6: At this point, it is worthwhile to revisit the 
assumption of 1 bit/s/Hz fixed transmission rate. According 
to the scheduling scheme, receivers select their transmitting 
nodes by feeding back their indices. Accordingly, the nodes 
transmit independently at 1 bit/s/Hz. Receivers decode their 
scheduled transmitters independently, by treating concurrent 
interference as noise. In other words, it is assumed that 1) 
the transmitters do not adapt their transmission rate to the 
instantaneous channel realizations; and 2) the receivers do not 
attempt to perform any interference cancelation. It is reason- 
able to expect that a higher throughput can be achieved if we 
allow rate adaptation and interference cancelation at the cost of 
more feedback overhead and higher computational complexity. 
However, what Lemma [3] tells us is that the return is at 
most the multiplicative factor log log n. Simulation results in 
Section [V] indicate a rather fast convergence to the asymptotic 
limits with increasing number of nodes. In a practical system 
with finite (but maybe large) n, the term log log n is a small 
number. On the other hand, given the decentralized scheduling 
policy adopted here, it is not straightforward to determine 
the adaptive transmission rate for each transmitting node (cf. 
Remark [TJ. 

IV. How Fast Can to Grow? 

As discussed in Remark [3j in Phase 1 there is a tradeoff 
between the number of relays to that serve as conduits between 
the source and destination nodes and the mutual interference 
caused by the transmissions. The same is true for Phase 2. 
This brings up the question: What is the optimal to that 
maximizes the throughput? This is equivalent to asking the 
maximum throughput of the network. In this section, we show 
that both hops of the proposed decentralized scheme succeed 
in achieving O(logn) throughput scaling (without taking the 
feedback overhead into account). We then quantify the feed- 
back overhead and conclude that, under the condition that the 
product of the block duration and the system bandwidth scales 
faster than log n log log n, the feedback overhead is negligible 
and therefore the useful throughput of the proposed scheme 
is given by 9 (log n). As a by-product in characterizing the 

4 If one schedules the transmission randomly, the average receiver SINR 
can be shown to be — —. r . 



throughput upper bound of the first hop, we also conclude 
that O(logn) is indeed the best throughput scaling even if 
centralized scheduling is allowed. Thus, as far as throughput 
scaling is concerned, operating the network in a decentralized 
fashion, with local CSI at the receivers and low-rate feedback 
to the transmitters, incurs no loss. 

A. Phase 1 

Earlier, in Section IIII-AI the lower bound (O of the system 
throughput of Phase 1 was found. This lower bound was 
adequate for the discussion in that section which assumed 
a large n and fixed to. However, in seeking to determine 
how the throughput scales with to, the lower bound (O 
might considerably underestimate the true throughput. In light 
of this, in order to address the question of optimal to, we 
reason as follows: First, we consider a genie-aided scheme 
by relaxing the assumptions of decentralized relay scheduling. 
Thus, the throughput scaling for a genie-aided network with in 
relays serves as an upper bound on the proposed decentralized 
scheme. We show that the throughput scaling law of this 
genie-aided scheme is 0(log?i). Next, we show that the lower 
bound (0 of the proposed decentralized scheme also achieves 
the O(logn) scaling. Thus, we are able to conclude that the 
throughput scaling of the original scheme of Phase 1 is given 
by O(logn), and is optimal in a scaling law sense. 

1) Phase 1: Upper bound due to genie-aided scheme: 
In this subsubsection, we establish the upper bound on the 
throughput scaling of Phase 1 based on the following genie- 
aided network. The genie-aided network has access to the 
full CSI of the network, and can coordinate the operation 
of the entire network, i.e., centralized scheduling is allowed. 
Therefore, the genie network can always achieve the maximum 
throughput in that, by assumption, it can enumerate all possible 
combinations of source-to-relay transmissions. Nevertheless, 
we still assume that independent encoding at the source nodes 
and independent decoding at the relay nodes. These constraints 
are needed to keep the genie-aided upper-bound result not 
too loose with respect to the proposed decentralized network 
(cf. Section |H}. Note that given a set of channel realizations, 
the successful source-relay pairs, in the proposed decen- 
tralized scheme, must also be successful in the genie-aided 
scheduling scheme. Thus, the throughput of the genie-aided 
scheduling scheme upper-bounds the proposed decentralized 
scheme. 

Theorem 2: Under the assumption of independent encoding 
at the source nodes and independent decoding at the relay 
nodes, one cannot achieve j§ £ §+2 throughput with probability 
approaching one. Conversely, with probability approaching 
one, (l-e)fk^ + 2 throughput is achievable for all e € (0,1). 

Outline of proof: The upper bound result is established 
by showing that, given m = + 2 relays, with probability 
approaching 1, one cannot find ^|^- + 2 sources whose con- 
current transmissions to the relays are all successful, even by 
enumerating all possibilities of choosing sources and mapping 
sources to relays. The achievability result follows from the 
fact that, by exhaustive search, with probability of one, one can 



find successful concurrent transmissions from (1 — e) 2 °f g " 2 + ^ 
sources to to = (1 — e) J"^ " 2 + 2 relays. 

See Appendix [B] for detailed proof. □ 

Remark 7: These results may be of interest in their own 
right, since, given the assumptions of independent encoding 
at the transmitters and independent decoding at the receivers 
(reasonable assumptions in ad hoc networks in which global 
CSI is not available to enable cooperative encoding and/or 
decoding), Theorem |2] establishes the upper bound and the 
achievable throughput scaling that are valid even if centralized 
scheduling is allowed. The 8(logn) throughput exemplifies 
the interference-limited nature of the network, and this sub- 
linear throughput scaling (compared to other works, e.g., [3]- 
[5], [13]) precisely demonstrates the price one has to pay for 
not having global CSI knowledge to mitigate the interference. 
Indeed, recent works show that if one allows either cooperative 
decoding between receivers (see, e.g., [5]) or cooperative 
encoding between transmitters (see, e.g., [13]), one can indeed 
avoid/cancel interference, enabling linear throughput scaling. 

2) Achievable throughput scaling of Phase 1: Theorem [2] 
states that, with high probability, there exists a valid group 
with to = (1 — e) 210^2 + 2 sources such that all transmissions 
are successful. However, the proof is nonconstructive: it does 
not afford insight into how to find such a set in practice. 
The proof assumes that there is a genie with global channel 
information that can enumerate all possibilities and select 
a good one for scheduling. In contrast to the genie-aided 
scheme, the opportunistic relaying scheme seeks to operate 
in a decentralized manner, and it is not clear whether this 
operational simplification incurs a loss in the scaling order of 
the throughput. Serendipitously, it can be shown that the log n 
scaling is also met by the lower bound in ((7). To see this, we 
examine the asymptotic behavior of (0. 

Consider the exemplary case of to = logn and s = logn — 
log log n. With n — > oo, the term "("- 1 )-^(»- m + 1 ) _> \ The 
term (l — (1 — e _s )™) is independent of to, and approaches 1 
for s = logn — log logn as 71^00. Therefore, a throughput 
of O (logn) can be achieved as long as Fy(s — 1/p) = O (1). 
Indeed, for m = logn, the interference term Y in ©, by 
the central limit theorem, can be approximated as Gaussian 
random variable with mean and variance both equal to log n. 
Now, we have 

F Y (log n - log log n - 1/p) w F Y (log n) = - , (14) 

due to the symmetry of the Gaussian distribution. Conse- 
quently Ri « i log n. This result implies that for m = log n 
relays, each running the two-hop opportunistic relaying pro- 
tocol, it is possible to schedule up to logn source nodes to 
transmit simultaneously, but half of them will fail to satisfy 
the SINR requirement due to the multiple access interference. 
In terms of throughput, this example yields \ log n, which 
confirms that the scheme is in fact order-optimal in achieving 
a throughput of Q (logn) at Phase 1. 

B. Phase 2 

In this subsection, we will show that the optimal value of 
to in Phase 2 exhibits a sharp phase transition phenomenon. 
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That is, m = log " lo ^ o 1 °| n 1 ^ PR + 1 succeeds in retaining the 

linearity of R 2 in m, but m = log " +lo f o 1 ° s 2 "~ 1/p " + 1 does 
not. As far as the scaling law is concerned, this implies that 
the throughput of Phase 2 scales as O(logn). 

Theorem 3: For Phase 2 of the two-hop opportunistic relay- 
ing scheme, if the number of relays m = lo s " ~ lQ g|°g " ~ 1 / p " _|_ 
1, then i?2 = (to) = 9 (logn). Conversely, if m = 

logn+loglogn-l/p fl + ^ ^ ^ = Q ^ m y 



log 2 



Proof: For convenience, we repeat R 2 of (fTTb and ( 112) : 

i? 2 =m(l- (F(l)) n ) (15) 

e -l/pH\" N 



m 1 - 1 



(16) 



log n — log log n— 1/pR , -i i_ 

With m = — + 1, we have 



l-F(l) 



log 2 



,-(m-l)lo g 2-l/ PJi 



logn 



Then, 



(F(l))"= (l-^V = 



e — log n+Q( ° B „ " ) _ g— logn+o(logn) 



(17) 



where we have used the fact that, for small x, log(l — x) = 
—x + 0(x 2 ) and e x = 1 + 0(x). Thus, most of the transmis- 
sions meet the SINR threshold (with probability 1 — 0(l/n)), 
and consequently the throughput R 2 is given by m(l — 0(—)). 
R 2 = 0(to) = O(logn) follows readily. 

Similarly, when m = l °s n+l °s^°s n - 1 / pr _|_ ^ we jj ave ^ 

F(l) 



n log n 



and 



(*■(!))" 



log 2 



= 1 - 0(1/ logn). 



(18) 



Now, in contrast to the case of m = log " lo ^" x l pR + 1, 

when we increase m to los " +lo f o '°g"~ 1/pfi + 1, Phase 2 of the 
two-hop scheme cannot support a throughput that scales with 
m. With probability one, the SINRs cannot meet the threshold. 
In this case, the throughput does not scale linearly with m 
anymore, i.e., R 2 = o(m). □ 
This G(log??) scaling result is consistent with the random 
beamforming scheme of [14], an outcome that is not surprising 
in light of the connection discussed in Remark [4] 

C. Feedback Overhead Analysis 

One of the contributions of the paper is the proposal of a 
two-hop scheme that alleviates the assumption of full CSI at 
the transmitters and the assumption of centralized scheduling. 
In the proposed scheme, only CSI at the receivers is employed, 
but low-rate feedback from the receivers to the transmitters is 
assumed to enable scheduling. In this subsection, we quantify 
the overhead due to feedback, and formalize a sufficient 
condition for which the overhead is negligible. 



1 ) Overhead per fading block: 
Feedback overhead in the first hop: In the first hop, each relay 
schedules one source. Since a total of n sources need to be 
identified, the feedback overhead per relay is log 2 n bits. The 
overall feedback overhead of all relays is thus given by m ■ 
log 2 n, which scales as 0((logn) 2 ) since m = 9(logn). 

Feedback overhead in the second hop: In the second hop, any 
destination feeds back the index of a relay only if there is one 
relay meeting SINR > 1; otherwise, no feedback is sent. One 
needs 0(log 2 m) = 9 (log log n) bits to identify a relay. The 
number of users that capture a good SINR, and consequently 
feed back follows the binomial distribution Bi(n, q) with 
q being the probability that the destination will provide a 
feedback. Then, the average overall feedback overhead is given 
by the average number of destinations that feed back, nq, times 
the number of bits of each feedback. 
To calculate q, we have, 

q = Pr[the destination feeds back index] 
= Pr[U™ =1 {the rth relay has SINR > 1}] 

m 

< ^2 Pr [ the rth rela y has SINR > 1] 

= m(l-F(l)) 

= ^re-Vp, (19) 



where the last equality is due to ([Toi l. The quantity in (fT9l is of 
the order of 9 ((log nf /n) when m = log "~ lo f o 1 ° s 2 "' 1/p " + 1 
(cf. Theorem 0. 

Finally, the average feedback overhead of the second hop 

is O(ngloglogn) = O ((log n) 2 log log n). 

2) Condition for 9(logn) useful throughput: In this work, 
a quasi-static fading model is assumed. Specifically, it is 
assumed that channel gains are fixed during the transmission 
of each hop, and take on independent values at different hops. 
The feedback overhead per block is analyzed above. The 
throughput of the system scales with the duration of the block 
as well as system bandwidth. That is, the total throughput 
scales as @(TW ■ logn). 

Now, we can find conditions of TW for which the feedback 
overhead does not imperil the throughput in the sense of 
throughput scaling. This is true when the feedback overhead 
is less than (in terms of scaling) the total throughput of each 
hop. Let us look at the second hop, which has larger feedback 
overhead. To have 9 (logn) useful throughput, we need 

(logn) 2 log logn = 0(TW ■ logn), 

which holds whenever TW = f2 (logn log logn). 

In practical system design, T can be as large as the 
coherence time of the channel T c , and W can be as large 
as the coherence bandwidth of the channel W c . For a typical 
wireless channel, this condition is easy to meet. This is due 
to the fact that typical wireless channels are underspread, that 
is, they satisfy T C W C S> 1. In typical urban environments, 
the coherence bandwidth is of the order of several MHz, and 
the coherence time is of the order of milliseconds [16]. Thus, 
the product of coherence bandwidth and coherence time is 
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Fig. 2. First hop average throughput Ri as a function of the number of relays for n = 1200 S-D pairs. From the top: simulation results utilizing all source node 
assignments, simulated results with distinct scheduled nodes, lower bound with optimized threshold s, and lower bound with threshold s = logn — log log n. 



of the order of 10 3 . As a concrete example, suppose the 
carrier frequency is f c = 900 MHz, and the delay spread is 
T c i = 1 /is. Based on the definitions of coherence bandwidth 
and coherence time in [16], the coherence bandwidth is given 
by W c = l/{2T d ) = 0.5 MHz. The coherence time depends 
on velocity v, where let us assume v = 3 km/h. This leads 
to a maximum Doppler spread of D s = f c v/c = 2.5 Hz, 
and accordingly, to a coherence time of T c — 1/AD S = 100 
ms. In this example, T C W C = 5 x 10 4 , which makes T C W C > 
log n log log n hold even for extremely large n. For example, 
for n = 1.0 x 10 s , logn log logn = 53.7. 

D. Two-Hop Communications 

With the results in previous subsections, we can conclude 
the achievable throughput scaling of the scheme in the follow- 
ing theorem. 

Theorem 4: Under the setup of Section HU and given 
TW = £!(log n log log n), the proposed two-hop opportunistic 
relaying scheme yields a maximum achievable throughput of 
6 (logn). 

Remark 8: The throughput scaling results in this paper 
afford a multiuser diversity interpretation. To see this, it is 
useful to take a closer look at the first hop. The power of the 
signal of each scheduled link is given by logn [8] (due to 
multiuser diversity). In the regime where m is fixed (cf. Sec- 
tion UHJ, the signal power can mitigate the interference power 
(which is of the order of one), and therefore each scheduled 
transmission is successful. This translates to m bits/s/Hz total 
throughput of the first hop, as shown in Corollary Q] In the 
limiting operating regime with m = 0(logn) relays, the 
aggregate interference for each scheduled link is of the order 
of 0(logn). Now, the network saturates as the interference 
power is of the same order as the signal power. The system 



throughput is calculated as 0(logn) -0(1) = 9 (logn), 
since one has 0(logn) concurrent transmissions, and each of 
them has successful probability of 0(1). Further increasing m 
results in a decreased probability of successful transmission. 
Referring back to Remark[5] we see that the optimal m (in the 
sense of maximizing system throughput) is given by the order 
of multiuser diversity. The interpretation of the throughput 
scaling in terms of multiuser diversity is discussed in more 
detail in [19]. 

V. Numerical Results 

In this section, we provide some numerical examples pro- 
duced by simulations of the proposed opportunistic relaying 
scheme under Rayleigh fading. Throughout these examples, 
the SNR for both hops is set at 10 dB (p = p R = 10 dB). 

We examine in Fig. |2]the average throughput R\ of the first 
hop of the protocol and its various lower bounds. The figure 
contains four curves. The two simulation curves were obtained 
by averaging throughputs over 2, 000 channel realizations. The 
"simulated curve was obtained using all assignments of 
source nodes, while the curve marked "simulated R\ with 
distinct nodes" represents only assignments of distinct source 
nodes. The other two lower bounds shown are computed with 
d2): one is obtained by optimizing (|7]i over s (numerically); 
the other lower bound is for s — logn — log logn. Three 
observations are noteworthy relative to Fig. |2] First, both 
the simulated throughput and the analytical lower bound (0 
exhibit linearity with respect to m, consistent with the analysis 
of Section ITlI-AI Second, it is observed that when m exceeds a 
certain value (in this case, 6), the throughput R\ starts to fall 
off. Noting that log 1200 = 7, this effect is consistent with 
the analysis in Section IIV-AI that established that the linear 
increase in throughput with the number of relays holds only 
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Fig. 3. First hop average throughput Ri, second hop average throughput R2, and average system throughput R as a function of the number of relays m 
for n = 1200 S-D pairs. 



as long as to is of the order logn. Third, the lower bound 
of i?i of (0 becomes loose when m grows. The development 
leading to (O suggests two possible reasons for this behavior. 
The first is that the computation of Pr [iV m ] is based on only 
distinct source nodes. However, the close match between the 
two simulation curves in Fig. [2] eliminates this possibility. 
It follows then that the bound is loosened due to the series 
of lower-boundings of Pr [ YJJ^Y ] l eacrm g to (O being too 
conservative. 

In Fig. [3] we illustrate throughputs i? x and R2, as well as 
the corresponding system throughput of the full scheme given 
by R = i min{i?i, R?}- As discussed in Section[II] the trans- 
missions over the second hop are destined to be successful, 
since they are scheduled based on SINR measurements at the 
destination nodes, whereas the transmissions over the first hop 
are not guaranteed to be successful since they are based only 
on SNR measurements. As a consequence, we observe from 
Fig. [3] that R\ is lower than R2, and is the bottleneck to the 
system throughput, i.e., R = \R\. In addition, we observe 
that the optimal number of relays for Phase 2 is consistent 
with the analysis of Theorem [3] in Section lTlI-BI Nevertheless, 
both Ri and R 2 display the linearity in to as predicted by 
Corollaries Q] and [2] in Section [TTTJ 

The total throughput of the two-hop opportunistic relaying 
scheme is shown in Fig.[4]as a function of the number of nodes 
n. We observe that the throughput exhibits the logn trend, 
as predicted by Theorem [4] In fact, the system throughput 
curve can be perfectly approximated by 0.36 logn. Since 
the system throughput is always limited by Phase 1, i.e., 
R = i min{i?i, R2] = we also plot two bounds of ^Ri 
for reference. More specifically, the genie bound + 1 

(cf. Theorem [2]) serves as an upper bound, and the ^ log n 
curve from ( fl~4b serves as a lower bound for the system 



throughput. In Fig. [5] the optimal value of the number of 
relays, to, is shown versus the number of nodes, n. Comparing 
the values of to from the curve, with the value '? s " +2, which 

' 2 log 2 ' 

is the bound on the number of relays for the genie scheme in 
Theorem [2] we observe that the optimal to is very close to 
that of the genie-bound. This explains why the scheme can 
harness large portions of throughput as promised by Theorem 
4. 

VI. Other Considerations: Relay Cooperation and 

Delay 

One of the key contributions of this work is to propose 
an opportunistic relaying scheme that features decentralized 
relay operations and practical CSI assumptions. In this section, 
we discuss the case in which relays are allowed to cooperate 
in encoding/decoding. In order to isolate the impact of the 
relay cooperation on the two-hop scheme, we leave the CSI 
assumptions unchanged. Specifically, it is assumed that the 
relays have full CSI knowledge of the source-relay link, but 
have only partial index-valued CSI knowledge of the relay- 
destination link via feedback. This discussion will help identify 
the fundamental limits of the opportunistic relaying scheme. 
Finally, we briefly address the issue of network delay. 

A. Cooperative Relays 

In the proposed opportunistic relaying scheme, we assume 
the relays perform independent decoding (in the first hop) and 
independent encoding (at the second hop). In particular, the 
relays treat the received interference as noise, and no attempt 
is made to cancel the mutual interference caused by concurrent 
transmissions. As a consequence, the system is interference 
limited. In this subsection, we address the question: How 
does cooperation between relays in decoding/encoding change 
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the scheduling operation and the throughput scaling? For 
example, it is conceivable that the relays could be implemented 
as infrastructure nodes that are connected to a wired backbone. 
This setup has been referred to as a hybrid network (see for 
example [20] and references therein). 

When the relays are allowed to fully cooperate in decod- 
ing/encoding, they can be considered to be a multi-antenna 
array. Accordingly, the first and second hops are equivalent 
to a MIMO MAC with receiver CSI, and a MIMO BC with 
partial transmitter CSI, respectively. Now, the scheduling in 
Phase 1 can be simplified. It is well-known that for the MIMO 



MAC, the sum-capacity can be achieved by allowing all 
users to transmit. The receiver can retrieve the data via some 
sophisticated signal processing algorithm, e.g., MMSE-SIC 
(minimum-mean square estimator with successive interference 
cancelation) [21]. The optimal scaling in the large n and fixed 
m regime is given by m log log n. However, if we seek to 
achieve only linear scaling in m, it suffices to schedule any 
m source nodes for transmission. With high probability, the 
resulting m x m channel is well-conditioned, and a spatial 
multiplexing gain of m is achieved [16]. In contrast, Phase 2 
does not benefit from the cooperation of relays. This is be- 
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TABLE I 

Dependence of Throughput Scaling on Relay Cooperation in Encoding/Decoding and CSI Knowledge of the Relay-Destination 

Link 



Scenario^ 


Throughput Scaling of Ri 


Throughput Scaling of R2 


Throughput Scaling of R 


Case 1 


0(logn) 


O(logn) 


0(logn) 


Case 2 


6(n) 


O(logn) 


e(logn) 


Case 3 


0(n) 


8(n) 


O(n) 



a Case 1: independent decoding/encoding at the relays; perfect CSI in the first hop and partial CSI in the second hop 
Case 2: cooperative decoding/encoding at the relays; perfect CSI in the first hop and partial CSI in the second hop 
Case 3: cooperative decoding/encoding at the relays; perfect CSI in the first hop and perfect CSI in the second hop 



cause, with only partial transmitter CSI, and since destination 
nodes are not allowed to collaborate, arbitrarily selecting m 
destination nodes cannot yield a throughput linear in m [14]. 

The impact of relay cooperation on throughput scaling 
exhibits similar behavior to that demonstrated for scheduling. 
Phase 1 benefits from relay cooperation, and in principle 
i?i = 8(n) is possible (note that the capacity of anxn MIMO 
channel scales linearly with n [22]); Phase 2 is still bounded by 
0(log7i), and becomes the bottleneck of the two-hop scheme. 
The reason for this is that, due to the lack of full CSI at the 
relays, there is no way to generate more than O(logn) parallel 
channels. The reader is referred to [14] for a discussion of the 
impact of CSI knowledge on MIMO downlink channels. 

We summarize the discussion of relay cooperation in Ta- 
ble U In the first two scenarios, we examine the impact of 
cooperation in decoding/encoding at the relays on system 
throughput by fixing the CSI assumptions to perfect receiver 
CSI in the first hop and partial CSI in the second hop. 
Specifically, case 1 corresponds to the setup considered in 
Section [II] where the relays perform independent decoding 
(in the first hop) and independent encoding (in the second 
hop), and case 2 allows for cooperation among relays in 
decoding/encoding. In the comparison, we also include the 
optimistic scenario, case 3, where the relays are assumed to 
have full CSI knowledge of both the source-relay link, and the 
relay-destination link and the relays are allowed to cooperate. 
In this case, 0(rt) throughput is obtained in both hops, a result 
not surprising from the MIMO theory [22], [23]. From the 
table, one can readily identify that CSI plays a critical role 
in determining if linear throughput scaling is achievable. This 
observation justifies our study on throughput scaling based on 
the seemingly pessimistic, yet practical, assumptions on CSI 
knowledge in this paper. 

It is important to point out that in the above discussion, the 
focus is on CSI. In cases where perfect CSI is available at 
the relays, but cooperative decoding/encoding is not available 
(e.g., due to nodes located randomly), different conclusions 
can be drawn. For example, one can operate the two-hop 
amplify-and-forward scheme [3] to achieve 0(n 1 / 2 ) through- 
put scaling. A detailed discussion of the case of perfect CSI 
but no cooperation between relays is outside the scope of this 
paper. It is also important to point out that the discussion 
applies only to the underlying Rayleigh fading model. For 
other fading models, the opportunistic relaying scheme may 
exhibit a different scaling law. See [24] for discussions of 
throughput scaling under more general fading models. 



B. Delay Considerations 

There is always a tension between opportunistic scheduling 
and delay considerations [8]. The delay issue is more salient 
in the two-hop scheme than in the cellular setup [8], because 
packets transmitted by one particular source in Phase 1 must 
be buffered at a relay, until that relay schedules the original 
destination during Phase 2. While one can partially relieve 
the problem by, say, prioritizing the destination in cases when 
a relay receives multiple requests from multiple destinations 
(including the destination of interest, of course), the delay 
may still be large. The detailed study of end-to-end delay is 
currently underway. 

VII. Conclusion 

In this work, we have proposed an opportunistic relaying 
scheme that alleviates the demanding assumptions of central 
scheduling and CSI at transmitters. The scheme entails a two- 
hop communication protocol, in which sources can commu- 
nicate with destinations only through half-duplex relays. The 
key idea is to schedule at each hop only a subset of nodes that 
can benefit from multiuser diversity. To select the source and 
destination nodes for each hop, relays operate independently 
with receiver CSI only, and with an index-valued feedback to 
the transmitter. The system throughput has been characterized 
for the operating regime in which n is large and m is relatively 
small. In this case, the proposed scheme achieves a system 
throughput of m/2 bits/s/Hz, while the upper bound with 
full cooperation among relays and full CSI is (m/2) log log n. 
Moreover, we have further shown that, given that the product 
of the block duration and the system bandwidth scales as 
fi(lognloglogn), the achievable throughput scaling of the 
proposed decentralized scheme is given by O(logn), which 
is the optimal scaling even if centralized scheduling is al- 
lowed. Thus, operating the network in a decentralized fashion, 
with only CSI at the receivers and low-rate feedback to the 
transmitters, incurs no penalty. Finally, compared to the linear 
throughput scaling results reported in the literature (see, e.g., 
[5] and [13]) with more optimistic CSI assumptions, this work 
quantifies the price that one has to pay for not being able 
to mitigate interference. The delay behavior of the proposed 
opportunistic relaying scheme is left for future work. 

Appendix A 
Characterization of Interference Y of © 

In this appendix, we characterize the statistical properties of 
the interference term Y in (0|. More specifically, it is shown 
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that, asymptotically in n, each individual term that comprises 
Y has an exponential distribution, and all interferers are 
asymptotically independent. It is also illustrated by numerical 
results that these asymptotic trends are achieved quickly, 
enabling the approximation of Y as a chi-square random 
variable with 2(m — 1) degrees of freedom. 

For notional convenience, we denote the channel connec- 
tions from n sources to the relay as X\,... ,X n . According 
to the scheduling of Phase 1, for each time-slot, i.e., each 
realization of Xi,. .. ,X n , the desired signal strength is the 
maximum among all connections. The interference term Y is 
the summation of (m— 1) out of the remaining (n— 1) channel 
connections. 

We first show that each interferer is asymptotically expo- 
nentially distributed. By the law of total probability, for events 
B and A, we have 

Pr[B] = Pt[B\A] Pt[A] + Pt[B\A] PrfA], (20) 

where A denotes the complement of the event A. Now define 
the event B as {X; < a;,}, and A as {Xi is not the maximum}. 
Then we have for the cdfs 

F Xi (x i )=F Xi (x i \A)Pi[A]+F Xi (x i \A)PT$\. (21) 

In our i.i.d. model, by symmetry, each node has probability 
of 1/n to be the maximum, i.e., Pt[A] = 1 — 1/n. Thus, the 
above equation can be written as 

FxM) = FxM\A) (1 - -) + F Xi ( Xi \A)- . (22) 

n j n 

Therefore, we have 

F Xl (xi\A) -> F Xi (xi ) as n -f oo, (23) 



and thus, asymptotically, each interferer is still exponentially 
distributed. 

While the above results are of an asymptotic nature, 
numerical result shows that they hold for practical values 
of n as well. For example, in Fig. |6j the empirical pdf 
f(xi\Xi is not the maximum) is plotted together with the pdf 
of Exp(l), i.e., f(x) = e~ x , x > 0, for various values of n. 
It is seen that the empirical pdf is well approximated by the 
standard exponential distribution. 

Next, we show that the interferers are asymptotically 
independent. Define A as {none of Xi,...,X m -i is 
the maximum}. The event A is then {at least one of 
X%, . .. ,X m _i is the maximum}. Again, by the law of total 
probability, 

F( Xl ,...,x m -x) =F{x 1 ,...,x m -i\A)Pr[A] 

+ F(x 1 ,...,x m - 1 \A)Pr{A\. 

Due to the underlying i.i.d. assumption, Pr[A] = (1— l/n)(l— 
1/(77,-1)) • ■ ■ (l-l/(ra-m + 2)) -> 1 and Pr[A] -> when 
n is large and m small relative to n. Then it follows readily 
from d24l i that in the regime of interest 

m— 1 

f(xi, . . . ,x m -i\A) -> f(xi, . . . ,X m -l) = J| f(Xi). (25) 

i=l 

Therefore, asymptotically, all interferers are independent. 

Combining the facts that, in the regime of interest, 1) each 
interferer is exponentially distributed and 2) all interferers 
are independent, the aggregate interference Y can thus be 
modelled as chi-square with 2 (to — 1) degrees of freedom. 
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Numerical result shows that this approximation is accurate for 
values as low as n = 40. 

Appendix B 
Proof of Theorem 2 



The linearity of the expectation yields 



Here, we prove Theorem [2] of Section IIV-AI For conve- 
nience, the theorem is repeated below. 

Theorem 2: Under the assumption of independent encoding 
at the source nodes and independent decoding at the relay 
nodes, one cannot achieve throughput with probability 

approaching one. Conversely, with probability approaching 
one, 5^3+2 throughput is achievable for all e € (0, 1). 

Proof: The proof relies on the probabilistic method [25]. 
The basic idea of the probabilistic method is that in order 
to prove the existence of a structure with certain properties, 
one defines an appropriate probability space of structures and 
then shows that the desired properties hold in this space with 
positive probability. This method of proof has been seen in 
various subjects of information theory, for instance, see [26, 
Ch. 8], which studies the bandwidth scaling problem in the 
context of spectrum sharing. The line of our proof follows 
[26]. 

The upper bound is established by the genie-aided schedul- 
ing, which performs an exhaustive search for the maximum 
concurrent successful transmissions. Specifically, in testing 
whether m bits/s/Hz is achievable, the genie-aided scheme 
enumerates all m-element subset of source nodes and tests 
whether the resulting m transmissions to the to relays are all 
successful. According to the genie-aided scheme, we define 
the probability space ft = {(A, 7r) : A C {1, . . . , n}, \ A\ = 
to, 7r is any permutation on {1, . . . ,to}}, where A denotes 
a random m-set of all n source nodes and ir denotes any 
possible TO-to-TO mappings from m source nodes in A to to 
relays. Let B\ be the event {all nodes in A can transmit 
simultaneously and successfully under mapping rule ir} and 
I\ the corresponding indicator random variable, i.e., I\ = 
if ,, > i Vi € A), where the subscript 

Rn(i) denotes the corresponding relay for source i under 
mapping rule tt. For any ir, we have {-/?„• (i),Vi € .4} := 
TZ = {1, . . . , m}. Finally, define the number of valid sets that 
satisfy the SINR threshold as X(m) = J2a I A- 
Then 



E[I A ] = Pr[i%] 



Pr 



SINR 



pi 



> 1, Vi e A 



= Pr 



SINR 



pi 



> 1 



(Pm) m , 



(26) 

(27) 



where 
SINR 



pi 



follows from the fact that for i,j € A, i ^ j, 



and SINR 



pi 



are i.i.d. The termp m = 1—F(i) 



in (IZ/D is the probability that a transmission is successful when 
there are to concurrent transmissions, and F(-) is the cdf of 
the SINR computed in ( ITOt . 



E[X(m)] 



to! (p n 



(28) 



Then, the upper bound is established by showing 
Pr[X(m) > 1] -> when m = + 2. This can be seen 
from Markov's inequality: 



Pr[X(m) > 1] < E[X(to)] 



(n 



{Pr, 



-i/p , 



< (np m ) m < (^^) m 

— m(logn-(m-l) log2-l/p) 

< e m(logn-(m-l)log2) 



(29) 



Now substituting to = H 

Pr[X(m) > 1] < 



2 into 



we have 



3 — logn+o(logn) 

= 0(i). (30) 

What ((30l) tells us is that when to = + 2, the proba- 
bility of finding a set of m nodes for concurrent successful 
transmissions decreases to zero as n increases. Since the 
transmission rate is fixed at 1 bit/s/Hz, it is equivalent to 
concluding that, with probability approaching one, + 2 
bits/s/Hz throughput is not achievable. 

Next we look at achievability. In proving the achievability 
result, we consider a variant of the genie-aided scheme used 
above. Here, the scheme divides the total n sources into to 
groups Qi,i = 1, . . . , to with each group having n/m sources. 
Each group is associated with one relay node. For example, 
as illustrated in Fig. [7] without loss of generality, we can 
label the total source nodes from 1 to n and assign sources 
{1, . . . , n/m} to Qi, {n/m + 1, . . . , 2n/m} to Q2, and so 
on. In testing whether to concurrent successful transmissions 
are possible, each relay chooses one source from its own 
group 2 Following the scheme, we define the sample space 
ST = {A : \A\ = m,R(i) ^ R{j) Vi,j 6 A), where R(i) 
denotes the index of the relay associated with the group to 
which source i belongs. Also define 1^ as the indicator random 
variable of the event {transmission from source i to relay R(i) 
is successful, Vz £ A}. Finally, let X'(m) = ^Aen' -^-4- 

To prove the achievability, we seek to find a lower bound 
on Pr[X'(m) > 1], or equivalently, an upper bound on 
Pr[X'(m) = 0]. We need the following probabilistic tool from 
[26]. 

Lemma 4: Let fj, = K[X'(m)] and A = 

T,Aen> E A' eat WaU']- Then, 

AnA'^o 



Pr[X(m) = 0] < e 



(31) 



Proof: Following the explanation in the proof of [26, 
Th. 10], the proof of the lemma follows in a fairly straight- 
forward way from [26, Lemma 7]. We skip the details for the 
sake of saving space. □ 

5 Similar scheme has been considered by Etkin [26, Ch. 8] in the context of 
characterizing the bandwidth scaling of spectrum sharing systems. Our setup 
is different from Etkin's scheme in that the number of nodes in each group 
is a function of m, which is not the case in [26]. 
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of IaIa' conditioned on {A G ft', A' G ft', \A n A'\ = q}. 
In the expressions (shown at the bottom of the page), d35| > 
upper-bounds (|34l by neglecting the interference coming from 
the sources belonging to A' H A in the second term of 
the product. In so doing, the two products Ilfee^i ■'■(■) an< ^ 
Y[ieA'\A^-(') m *E3 involve independent random variables 
now and therefore are independent (note that this is not true 
in (|34li). A minimal example is illustrated in Fig. [8] The upper- 
bounding in ( 1351 ) can be thought of as reducing the number 
of concurrent transmissions from m to m — q by keeping the 
elements {t : t G A' fl A} silent. The probability of successful 

transmission when there are to — q concurrent transmissions, 

-i/p 

denoted as p m -q, can be shown to be p m -q = 9 m- g -i ■ 



Fig. 7. Illustration of the genie-aided scheme used in the achievability proof. 

Source nodes are divided into m groups Qi,i = 1, , m, each with n/m 

nodes. Each group associates with one relay node and the nodes in the group 
have common receiver (i.e., the associated relay). A is formed by selecting 
one node from each group. 
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The proof of achievability is more involved than that of 
the upper bound. This is because /4's are generally not 
independent. In the upper bound proof, however, dependence 
among i^'s is irrelevant due to the linearity of the expectation. 
The quantity A in Lemma [4] is a measure of the pairwise 
dependence between the /4's. Note that, in the case when 
7/l's are all independent, the lemma reduces to Pr[X'(m) = 
0] = e _At , a result which can be reached by direct probability 
calculations. 

We begin with fx: 



/' 



ELY' (to)] = 1 

.)" 



''[^Aen' 1 ^] 



(-) (Pn 

\mJ 



(32) 
(33) 



where d32b follows from linearity of expectation, and d33l 
is due to the fact that, for any A G ft', all relays see 



i.i.d. channel realizations. The term p. n 



(cf. 423) 



denotes the probability of successful decoding when there are 
to concurrent transmissions in total. 

To compute A, let us start with computing the expectation 



Fig. 8. Example of A = {1, i} and A' = {1, j}, i 7^ j. We can upper-bound 
the SINR of source j in A 1 as YT^+^l 2 — T^p' wn i £n now ' s independent 
of the SINRs of source nodes in A. 

Now, we proceed with A. In particular, we have 



A = E E wm 



.Aen' ,4'efi' 

rn 

EE E * 

9=1 Aen' a' en' 
\AnA'\=q 



A£Q' 

A' en' 
\AnA'\=q. 



= E - 

^— ' Vto 

q=l 



1 

q I Vto 



(p m ) m (p m - 9 r 



In order to apply Lemma [4] we check 

a _^ (?)(^-i) m " g ( Pm -,r-' _f^ 

u (*r (p^r hi 



where a„ 



(p m -,r 



7^ 



Now on defining & 9 



E 



Ago' 
.4' en' 
|.AnA'|=g. 



< 



< 



= E 



keA 
fce.4 
.fce.4 
.fce.4 



7fc,it(/c) 


<T 2 /-P 


+ SteA 7t,fl(fc) 




lk,R{k) 


C7 2 /-P 


+ J2teA lt,R(k) 








lk,R(k) 


CT 2 /P 


+ J2teA lt,R(k) 








lk,R{k) 


CT 2 /P 


+ YlteA lt,R(k) 




t^k 


\m- 


Q 



>1 ' 



> 1 



>i) n 1 

>eA'\A 

>i) n 1 



li,R(i) 



a 2 /P + J2teA' lt,R(t) 



> 1 



l£,R(£) 



a 2 jP + J2teA'\A lt,R(l) 



> l 



> l 



E 



n 1 

UeA'\A 



ll,R{l) 



<? 2 /P + EteA'\Alt,R(e) 



> 1 



(34) 



(35) 



(36) 
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a q+ i/a q , we have that b q = ^szrrj^^T) 2 2 ^ T " 9 1 ' decreases 
with g, 6 g < 61. Therefore, a g < b\~ a\. 

On setting m = (1 - e)^j^ + 2 for any e e (0, 1), we 
have 

61 - 2(^-1) 2 

_ e log(m-l)-log(^;-l)+2(m-2)log2+l/p-log2 
- p-elog"+o(log™) 



1 



Furthermore, with m = (1 — e) 



A 
/' 



log n 
2 log 2 

1 ^ a l 



2, we have 



2=E a ^ ai E 5 ? 1 - y^t" 

q=l q=l 

(TK^- 1 )" 1 " 1 (Pm-i)""' 1 

m 2 (p™-!)" 1 - 1 1 



< 



n (p m ) m 1 - 61 

m 2 e l/p 2 2(m-l) 

n 1 — b\ 

p 2 log m-log n+2(m-X) log 2-log(l-f>i) + l/p 
-elogn+o(logn) 



o 



(1 



Finally, Lemma [4] yields 

Pr[X(m) = 0] < er n ' 
In words, ( f37l > tells us is that when to = 



(1- 



log n 



(37) 
2 for 



' 2 log 2 

any e 6 (0, 1), the probability of not finding a set of m nodes 
for concurrent successful transmissions decreases to zero as 
n increases. In other words, the probability of finding to 
concurrent successful transmissions with m = (1 — e) 



log 71 

2 log 2 



approaches 1, Again, since the transmission rate is fixed at 1 
bit/s/Hz, it is equivalent to concluding that, with probability 



approaching one, m = (1 — e 



log n 
2 log 2 



2 bits/s/Hz throughput 



is achievable by the genie-aided scheme. 

This achievability result, together with the upper bound ( f30b . 
completes the proof of the theorem. □ 
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